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The motion of massless spinning test particles is investigated using the Newman- Penrose 
formalism within the Mathisson-Papapetrou model extended to massless particles by 
Mashhoon and supplemented by the Pirani condition. When the "multipole reduction 
world line" lies along a principal null direction of an algebraically special vacuum space- 
time, the equations of motion can be explicitly integrated. Examples are given for some 
familiar spacetimes of this type in the interest of shedding some light on the consequences 
of this model. 
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1. Introduction 

The motion of a massive spinning test particle in general relativity is described by 
the Mathisson-Papapetrou model"'^'^ consisting of 10 equations 
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r)Qal3 

^ =[PAur0, 

aru 

for 13 unknowns: the generalized momentum of the particle P" (4), the antisym- 
metric spin tensor S"'^ (6) and the unit tangent vector U" = dx^/dr;/ to the world 
line used for the multipole reduction (3). 

For the model to be complete, it must supplemented by three additional equa- 
tions which in the literature arc usually called "supplementary conditions" (SCs) . 
There are three standard choices for the supplementary conditions: 

1. Papapetrou-Corinaldcsi:^ S"^up = 0, where u is some arbitrary timelike vector 
chosen for convenience of calculation; 

2. Pirani:^ S°'I^U0 = 0; 

3. Tulczyjew:^ S'^^Pp = 0. 

Each of these is an algebraic condition restricting the spin tensor to have 3 fewer 
independent components. If the vector in terms of which the supplementary con- 
dition is expressed is timelike, this just means that the spin tensor is equivalent to 
a spin vector lying in the local rest space associated with this timelike direction. 
In a sense, only the Pirani and Tulczyjew SCs are "intrinsic" relations between the 
various unknowns of the model and have more credibility as physical conditions 
rather than just a mathematical convenience which is "coordinate dependent", as 
is the Papapetrou-Corinaldes SC. 

The extension of the Mathisson-Papapetrou model to the case of a null multi- 
pole reduction world line with tangent I has been considered by Mashhoon.® The 
equations of the extended model have exactly the same form as (1) and (2) but 
with U replaced by I and tu (the proper time parameter along the timelike world 
line) replaced by A (a parameter along the null world line such that dx^/dX = Z") 

^R^PpJl'SP^ , (3) 

= [PAl]"'^. (4) 

Eqs. (3) and (4) must then be solved assuming some SC. Consider only the more 
physically motivated SCs of Pirani and Tulczyjew, with Pirani's conditions now 
naturally generalized to S°'^ljs = 0. Furthermore, since we are interested in the 
massless limit of the Mathisson-Papapetrou equations, and since the mass of the 
particle is defined by m = — P • U, the massless limit would imply P • Z = 0, so we 
assume this to be true here. 

Mashhoon^ found Tulczyjew's condition S°'^P0 = to be inconsistent if P is 
also null: P • P = 0. While we could examine this condition for the nonnull case, 
it seems most reasonable to consider Pirani's SC as the only physically meaningful 
supplementary condition. Under these conditions Mashhoon^ has shown that the 
particle follows a null geodesic (so that one can assume A to be an afRne parameter) 
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with its spin vector (denoted by <S) 

5^ - \^'''"'^luSc,p (5) 

parallelly transported along the null world line with tangent I and either parallel or 
antiparallel to I. Moreover, again adopting Mashhoon's point of view, the general- 
ized momentum of the particle should be spacelike or even null and is determined 
by solving equations (3) and (4) supplemented by S^-^lfj = 0. The other compo- 
nents of the spin tensor not equivalent to the spin vector that one can introduce 
must also be determined.*^ Note that the spin vector depends on the scale of I and 
hence on the parametrization chosen for the null world lines to which it is tangent. 

In the present paper, using the Ncwman-Penrosc (NP) formalism, wc derive the 
equations of motion for a massless spinning test particle in any algebraically special 
vacuum spacetime for the special case in which the multipole reduction world line 
is aligned with a principal null direction of the spacetime. This leads to very simple 
equations which are straightforward to solve. Explicit solutions corresponding to 
some familiar Petrov type D and type N spacctimcs (including the Schwarzschild, 
Kerr, C metric, Kasner, and single exact gravitational wave spacetimcs) are derived 
and discussed. Finally, motion along null circular orbits is also studied and explicit 
solutions in black hole spacetimes are presented. 

We adopt the standard NP notation of Chandrasekhar^ which includes the met- 
ric signature H . 



2. Equations of motion in the Newman- Penrose formalism 

In this section we express the above equations for massless spinning test particles 
in the NP formalism. Let ei = l,e2 = n, 63 = m, 64 = m be an NP frame consisting 
of four null vector fields satisfying the relations I ■ n = 1 and m- fh = —1, so that 
the frame components of the metric tensor are 

fO 1 \ 

10 

-1 

VO -1 / 



(6) 



"It is worth noting that the Mathisson-Papapetrou model does not place a priori restrictions on 
the causal character of U and P and there is no agreement in the literature on this freedom. 
For example, Tod, de Felice and Calvani^ consider P to be timelike assuming that it represents 
the total energy-momentum of the particle while they do not impose any causality condition on 
the world line U, which plays the role of a mere mathematical "tool" to perform the multipole 
reduction; on the other hand according to Maslilioon,*' P can be considered to be analogous to 
a canonical momentum for the particle, in which case there is no reason to prefer a particular 
causality character for it, while the world line U must be timelike in the massive case or null in 
the massless case because it is a world line intimately associated with the test particle location. 
This uncertainty in the model itself, i.e., the uncertainty in the SC, has to be considered as a 
limitation of our description of spinning test particles in general relativity. We are aware of this 
problem and only consider the Pirani SC here. 
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Exactly like the NP representation of the electromagnetic Faraday 2-form, the 6 
real components of the spin tensor S are represented in terms of 3 complex scalar 
quantities Sq, Si, S2 



{Sab) = 



( 



\ 





'Si - Si 
-So 
-So 



Si + Si 


^2 
^2 



So 

-S2 



-Si + Si 



So 
-S2 
Si -Si 




(7) 



which have the well-defined spin- weights'' —1,0, 1 respectively.^ 

We first adapt the NP frame to the null world line of the multipole reduction 
by choosing its tangent as the first frame vector: I = ei. The zero mass condition 
= P ■ I = P"^ then forces the second component P^ of P to be zero (and hence 
P ■ P = — 2|P^p, so P is spacelike). Eqs. (3) and (4) then take the following form 
in a vacuum spacetime where the Riemann curvature tensor reduces to the Weyl 
curvature tensor (NP scalars V'o, V'i)V'2,V'3)V'4) 



dA 



TvSo — KS2 , 



2e& 



dPi 

"dA" 
dP3 

dA 



= -(e + e)pi - ttP^ - TtP^ - iPsSo - Va^o + 2^251 + 2^)251 - ^1^2 - Vi'52 



= fiP' + {e- e)P^ - 2^iSi + ^2So + 1P0S2 , = kP-" + kP 



(8) 



where D/dA = Z"Vq = D is the directional derivative along I. As usual the equation 
for is just the complex conjugate of the one for P^. 
Pirani's supplementary conditions are simply:'^ 



{Si + Si) = 0, So = 



(9) 



i.e., = is/2 is purely imaginary with s real and whose meaning is discussed below 
(don't confuse s with the symbol for the spin- weight of tensorial quantities). The 
only nonvanishing frame components of the spin tensor are then S24 = —S2 and 



'^An NP frame aligned with a given null direction / is only fixed up to boosts and rotations which 
leave the directions of I and n unchanged, i.e., up to class III tetrad Lorentz transformations 

A quantity C, which tranforms under this transformation as 

^ ^(P-l-9)/2gi(p-q)/2^ ^ 

is said to be a quantity of type (p, 5), with spin weight s = (p—(if)/2 and boost weight ui = {j>+q)/2. 
Most of the NP quantities have a well-defined spin and boost weight. The exceptions are the spin 
coefficients a, /3, 7, e and the directional derivative operators D, A, 6, 5 along the frame vectors. 
"^Mashhoon's notation used in® corresponds to / = S2 and g = s and his choice of frame transport 
law along I is equivalent to the conditions tt = and e -|- e = on the NP scalars which guarantee 
that I is the tangent to an aiBnely parametrized geodesic and that n is parallel transported along 
I. This latter assumption is unnecessary here. 
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= —is. With these conditions the equations of motion then become 

= Ks , = 12^5*2 , = iris - 2eS2 - , 

qA qA 

— = -(e + e)pi - 7rP3 - wP^ + is {1^2 - ^2) - ^'1^2 - V^i^s , 
qA 

= kP^ + kP^ , — = + (e - e)p3 + ^^^5 + ^^§2 • (10) 
dA 

The solution s = of the first equation would lead to 5*2 = and hence vanishing 
spin tensor, a trivial case that will not be considered here. One must therefore have 
instead k = 0, i.e., I is a null geodesic.^ Hence, these equations further reduce to 

^ = 0, ^=.is-2eS2-P^, 
dA dA 

dpi 

— = -(e + e)pi - 7rP3 _ ^p3 ^ _ _ ^^g^ _ ^^g^ ^ 
dp3 

— = (e-e)p3 + iV;iS + Vio^2 . (11) 

A further simplification can be made by using a type III rotation of the NP frame 
to set the spin coefficient e = 0. This corresponds to require that the congruence of 
I is affinely parametrized. The final set of equation is the following: 

_=0, —=n^s-P , 

— = -TtP^ - TtP^ + is (V'2 - i>2) - ^182 - i>lS2 , 

dp3 

— = iiJis + ^082 . (12) 
dA 

These equations are still general for any vacuum spacetime. In the special case that 
the 4-momentum is a null vector, equivalent to P^ = as seen above, then P = P^l 
and the previous set of equations specialize to 

ds _ d82 

_ = 0, ^=^^s, 

dPi 

= is (■i/'2 - ^"2) - i/'i'S'2 - ^1 "52 , = -iV'is + Vo'S'2 • (13) 
3. Interpreting the spin components 

The spin vector defined in Eq. (5) as the contracted dual of the spin tensor 

Sn^*gni'l^ _^ S = i{Si - Si)l + iSom ~ iSom , (14) 

when the Pirani supplementary conditions Si = is/2 and So = are imposed, 
reduces to the simple form 

S^ = sF, (15) 
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SO that s represents the heUcity of the particle (see p. 80 of reference 10 for the 

relation between the hclicity and spin- weight of a masslcss field) . Since s is constant 
along the world line, and I is parallelly transported along itself, the spin vector is 
parallel transported along recovering Mashhoon's result. 

However, the spin tensor is not completely determined by the spin vector alone 
as in the timelike case but also requires the independent component 52 

S = I A {S2m + S2fh) + ism Am . (16) 

The real vector {S2m + S^m) can be interpreted as a sort of "orbital angular 
momentum" contribution to the total spin (sec reference 10, p. 68), as measured 
with respect to the null world line of the multipole reduction associated with /. In 
a sense, since I is a null vector, the multipole moments taken with respect to any 
spacetime point on the null world line are always those of a particle in motion with 
respect to any physical observer. 

Alternatively, there are no observers such that the associated electric part of S 
vanishes identically as it does in the timelike case. To see this let us consider the 
Lorentz frame naturally associated with the NP frame, defined by 

I = -^{Eo +Ei), n= -^{Eo -El), m= -^{E2 + iE^) . (17) 

Using the notation U = Eq- Eq. (16) then becomes 

S = U AL{U) + *^^'>S{U) , (18) 

where 

L{U) = Re [521^2 - Im [S2]E3 , 

S{U) = Re [S2]E3 + Im [S2]E2 + s Ei , (19) 

which are respectively the electric and magnetic parts of S with respect to U. The 
notation for the spatial dual of a spatial vector X (spatial with respect to U, i.e., 
X-U = 0), 

l*(U)x]^l3 = , (20) 

as well as for the exterior product of two spatial vectors X and Y (spatial with 
respect to U) 

[X xr;y]'> = /"^'^[/^x^i; , (21) 

where rja^js is the unit volume four form, will be used here and below. Note that 
L{U) + iS{U) = V2S2 m + isEi . (22) 
The two algebraical invariants of S can then be written as 

h = Is-^S^f, = -L{Uf + S{Uf = s\ I2 = ls"^*S^0 = L{U) ■ S{U) = .(23) 
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The spin tensor can be classified in the same way as the electromagnetic 2-form. If 
the (complex) invariant 

I = + i*S''l'){S^0 + i*S^0) = h+ ih = 4(^0^2 - Si) (24) 

is nonzero, the spin tensor will be said to be non-null or nonsingular, while if it is 
zero the spin tensor will be said to be null or singular. In the present case since 
/2 = 0, it reduces to / = and the spin tensor is singular only when s = 0. Note 
that the spinning test particle model holds when a dimensionless quantity formed 
from s is very small: otherwise modification of the background due to backreaction 
should be considered and the particle cannot be considered to be a test particle. Due 
to the choice of P supplementary conditions, the quantity ^2 does not contribute 
to the spin invariant; therefore, requiring the magnitude of s to be small does not 
give any restriction on it. 

The 4- velocity u of any physical observer can be expressed in terms of its relative 
velocity with respect to U and vice versa^^ 

u = -t{U,u)[U + u{u,U)] , U = -i{U,u)[u + u{U,u)] , (25) 

where 

\HU,uW = \Hu,U)\\\ ^(U,u) = ll^l-\HU,uW = ^{u,U) . (26) 

Then 

S' = uAL(u) + *(")5(w) , (27) 

with L{u) and S{u) related to L{U) and S{U) by a boost transformation exactly 
like the electric and magnetic fields^^ 

P{U)L{u) = 7(C/, u) [L{U) + v{u, U)xu S{U)] , 

P{U)S{u) =j{U,u)[S{U) + u{u,U) XuL{U)] , (28) 

where P{U) = 1 + U ^ U projects orthogonally to U. Since Ii = 0, it follows that 
L{u) and S{u) are orthogonal, i.e., any observer will see the electric and magnetic 
parts of the spin tensor to be orthogonal. Moreover, from (19), it easy to see that 

L(U) + Eix S{U) = , (29) 

which compared with (28), gives L{u) = only for an observer moving with respect 
to U with velocity iy{u, U) = Ei, i.e., with unit speed, so that it does not correspond 
to any physical observer. 

4. Vacuum algebraically special spacetimes 

The Goldberg-Sachs theorem^^ states that a vacuum metric is algebraically special 
if and only if it contains a shearfree null geodesic congruence with tangent I (taken 
as the first vector of an NP frame), i.e., 

K = = a, -00 = = ■(/'I . (30) 
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Hence, assuming that the (single) line of multipole reduction associated with I is a 
world line in such a congruence, i.e., I is a principal null direction of the spacetime, 
then the set of equations (12) simplifies to 

_=0, -J^=^^s-P , 

H pi _ - d P3 

2_ = -^p3-^p3 + is{^j^-M , ^=0- (31) 

The formal integration of these final equations is straightforward, leading to 
s(A) = s(Ao) and P^{X) = P^(Ao) being constant along I and 

52(A) = is(Ao) / 7r(A')dA' - P\Xo){X - Ao) + &(Ao) , 



pi(A) = -p3(Ao) /" 7r(A')dA'-p3(Ao) /" 7f(A')dA' 

+is{Xo) I [MX') - ^2(A')] dA' + P^(Ao) . (32) 

JXa 

In the special case P'^(Ao) — (i.e., P null and aligned with I), the solution (32) 
reduces to 

52(A) = is{Xo) I 7r(A')dA' + S2(Ao) , 

P\X) = is{Xo) / [MX') - MX')] dA' + P^(Ao) . (33) 

Eqs. (32) and (33) suggest the presence of a "hehcity-connection coupHng" term 
in S2, as well as a "helicity-curvature coupling" term in P^. We will come back to 
these couplings when discussing explicit examples. 



5. Explicit examples in the type D case 

In each case we adopt a well-known NP frame with the above properties, assuming 
that the world line of the multipole reduction is an afSnely parametrized integral 
curve of the principal null direction I. Wo investigate the principal null direction 
corresponding to outgoing radial null geodesies in the Schwarzschild case and those 
which are asymptotically radial outgoing null geodesies in the Kerr case, and those 
in related spacetimes which are analogous to these. The ingoing null geodesies of 
this type can be treated similarly. 

5.1. The Schwarzschild case 

Consider the Schwarzschild solution^^ written in standard spherical coordinates 

ds^=(l-^)af-(l-^) ' dr^ - r\de^ + sin^ edcl>^) . (34) 
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Introduce the NP frame adapted to the outgoing radial null geodesies 



2M 
r 



dt + dr , n = 



1 



m = 



V2r 



sin 6 



'<t> ' 



(35) 



which has the following nonvanishing spin coefficients 



1 



2r 



2M 



M 



V2 



— — ) ' ^=^' /?=-a = ^cote,(36) 



4r 



and the single nonvanishing Weyl scalar ■02 = —M/r^. The integral curves of I can 
be expressed in the form 



t = r -rQ + 2M\n 



r~2M 



to , 



A - Ao + ro 



^0 , 



ho , (37) 



ro - 2M 

where to, ro, ^0,00 are arbitrary integration constants taken at A = Aq. Then 
Eqs. (10) reduce to 

dS2 dPi „ dP3 



dA=« 



dA 







dA 







(38) 



with solution s(A) = s(Ao), 52(A) = 5'2(Ao) - P^(Ao) (A - Ao), P\X) = P\Xo) and 
p3(A) = P^(Ao). The case P^{X) = is clearly trivial. 



5.2. The Taub-NUT case 

The Taub-NUT solution^^"^^ is often written in the form 



ds' = 



£2 



{dt + 2icosed(l)f -{r'^ +f) 



dr2 



+ d6»2 + sin^ edf 



,(39) 



where Atn = — 2Mr — i"^, and it represents the spacetime of a black hole with 
an additional gravitomagnetic monopole charge associated with the parameter £. 

Introduce the NP frame generalizing the previous Schwarzschild situation to 
which it reduces when £ = 



Atn 
1 



-dt + dr 



"=2 



dt 



\TN 



dr 



m = 



V2{r + i 

with nonvanishing spin coefficients 
1 1 



-2iicot0dt + d0 



sin^ 



(40) 



P 



(3 = -a = 



r — U ' ^ 
\f2 cot( 



iTN 



2{r^ii){T-itf 



i l + iM 
'2{r-i£f 



4 {r + i£) ' 



(41) 
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and the single nonvanishing Weyl scalar i/»2 = {£+iM) / {i+ir)^ . The integral curves 
of I have the form 

r^-2Mr-f 
t = r - tq + Mln- 



'0 



2Mro - £2 



-2 + £2 



-^^f ) ^( 

V\/M2+£2y V\/M2+£2 



r = A - Ao + ro , 6» = 6lo , <A = <Ao , (42) 

where to: '''O: and (po are arbitrary integration constants taken at A = Aq. 

The formal integration of (32) is then straightforward, leading to s(A) = s(Ao) 
and P^{X) = P^(Ao) being constant along I and 

^2(A) = -p3(Ao)(A-Ao) + 52(Ao) , 

olM^ _ . ^ r - ^0 + ^0 - - (^^' + (rp - M)2 - + ^2) 

P (A) - £.(Ao) I [(A-Ao + ro)2 + ^2p ^^^2 

+P\Xo) . (43) 

When Ao = ro = r+ = M + \/M2 + ^2 and -P^(Ao) = 0, i.e., the particle starts 
moving from the horizon r_|- , we have A = r and the solution is given by s(A) = s(Ao) 
and -P'^(A) = P^(Ao) constant along / and 

S,{r) = S,{r+) - P3(r+)(r - r+) , P\r) =ts{r+) j^^^ • (44) 

If in addition P''(r+) ~ 0, then S'2(r) = 5*2 (r+) and P is null and aligned with I 

P = ^^ir^)^^^.l- (45) 

The behavior of the varying component P^ of the total 4-momentum as a function 
of the radial parameter is shown in Fig. 1. Moving from r+ out to infinity the 
proportionality factor with I passes from zero to a maximum positive value at the 
real solution of — 3Mr2 — ir(^ + Mi"^ = and then returns asymptotically to 
zero. 



5.3. The C metric case 

The C metric solution, ^^"^^ which represents the spacetime associated with a uni- 
formly accelerating black hole, is often written in the form 

ds^ = Hdu^ + 2dwdr - 2Ar'^ sm^ Odrde - ^^^^J^d6'2 - r^Gd(t>'^ , (46) 



where 



H=l- — - A^r'^G - 2Acos6i[r - 3M(1 - ^Ir cos6')] , 
r 

G = sin^ e - 2AM cos^ 6 . (47) 
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-0.06 



Fig. 1. In the Schwarzschild, Taub-NUT and Kerr cases the varying component of the total 
4^ momentum is plotted versus the radial parameter r/r^, where is the (outer) horizon in each 
of these spacetimes. The initial conditions are chosen so that P^(r+) = with Ao = = r+ for 

all cases, s(r+) = 1 and e/M = 1 for Taub-NUT, and do = tt/2, Im[P''*(r+)] = l/{10-/2a) and 
a/M = .5 for Kerr. represents the component of the generalized moineutum of the particle 
along the world line of the multipole reduction which is assumed to coincide with a principal 
null direction of the spacetime. In the Scliwarzschild case this quantity is constant all the world 
line and identically zero for our choice of initial conditions. In the Taub-NUT case it varies but 
asymptotically returns to its initial value. In the Kerr case, it varies and does not return to its 
initial value, although it does approach the asymptotic value — l/(10r+), with our choice of initial 
conditions. 

We limit our considerations to the range of parameters M < La = 1/ (3^\/3) where 
the C metric is static. Introduce the NP frame 

I = -Or, n= , 

m = j^Or 1= -de + —= — , (48) 

^/2 V2rsin6' V2rVG 

with nonvanishing spin coefficients 

1 H AVG ^ V2cos0(l + SAM cos 9) 



/2 

[2ArG -h cos 6'(f -I- SAM cos 9)] , 

4rv G 

A^r^G + MA^r"^ cos^ BiA cos 6'r - 31 + M - A cos Or'^ 
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and Weyl scalar ijj2 = —M/r^. 

In order to write down the equations for the integral curves of Z, it is convenient 
to introduce the new variable 

Ar 



W = 



1 — Ar cos 6 



(50) 



in place of r [6 is constant along I), which has proven to be helpful in studying the 
horizons of the C metric, i.e., the solutions of the equation H = 0, equivalent to 



W^-W + 2MA = . 
In the case M < La its solutions are 

Wi = 2U , W2 = -U + VSV , Ws 

where 

M ' 



-u-Vsv 



U + iV = 



1 



V3 I La 



M2 



1/3 



I.e., 



y 1 /I M 

U = ;= cos -arccos-— 

V3 La 

In terms of these quantities one finds 



y 1 . /I M 

V = —= sin -arccos-— 
V3 \S La 



(51) 



(52) 



(53) 



(54) 



Wo + 



Wi 



(Wi-W2){Wi-W3) 



In 



W-Wi 



Wo-Wi 



W2 



(Wi~W2)iW2-W3) 

r = -X + Xo + ro, = Oq, 



In 



W -Wo 



W0-W2 



(W2-W3){Wi-W3) 



In 



1^0-1^3 



where 



Wo 



Aro 



Aro cos 6 



(55) 



(56) 



Finally, the formal integration of the equations of motion is straightforward, 
giving s(A) = s(Ao) and P^{X) = P^(Ao) constant along I and 



52(A) = i7Tos{\o){\ - Ao) - P'(Ao)(A - Ao) + S2{Xo) , 
P\X) = -27ro(A - Ao)Re[p3(Ao)] + P^Xq) , 

where in this case tt, depending only on 9, is constant along I. 
In the equatorial plane {0 = ■7r/2) we have ttq = —Aj^fl. Thus 

iA 



(57) 



52(A) 



2\/2 



s(Ao)(A - Ao) - P^(Ao)(A - Aq) + 52(Ao) 



P\X) = V2A{X - Ao)Re[p3(Ao)] + P^(Ao) . 



(58) 
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5.4. The Kerr case 

Consider the Kerr metric^^ expressed in standard Boyer-Lindquist coordinates 



ds^ = -{dt-asm^edc^f '''' 

5^ 



[(r^ + a2)d0 - adt]^ - — dr^ - Sd^^ ^^59) 
2-i ZA 



where the standard notation A = + — 2Mr, S 
adopted. Introduce the Kinnersley frame 



cos^ 9 has been 



m = —pz r— :- {ia sin^ Odt + sin ^9^1 + id^) . 



\f2 sin ^(r + ia cos 
The nonvanishing spin coefficients are 
1 1 



(60) 



ia sin Q 



a = — 

TT = 



r — iacost) ' 2 {r — ia cos 9)Ti ' \/2S 
i[Mi(r + iacos9) + a{rcos9 — ia)] 
2T,{r — ia cos 9) ' 
1 cos 6'(r + ia cos 6*) — 2m 



2\/2 sin 9 {r — ia cos 9)"^ 
ia sin 9 



cost 



y/2{r — ia COS 9)"^ ' 2-\/2sin0(r + iacos^) ' 

and the Weyl scalar is V2 = —M/{r — iacos9)^. The integral curves of I are 

J.2 _ 2Mr + 



(61) 



t = io + (r--ro)+Mln 

2M r 

, arctanh 

V'M2 - a2 L 

r = ro + A - Ao, 9 = 9o , 

a 



- 2Mro + 



r - M 



VM2 



arctanh 



V v'M2 - a2 
r - M 



arctanh 



/ ro-M 



VM2 



— arctanh 



M 



V\/M2 



(62) 



and represent outgoing null geodesies which are asymptotically radial far from the 
horizon. Eqs. (10) then have the following solution 

s(A) = s(Ao) , 
52(A) = 52(Ao)-P^(Ao)(A-Ao) 



—^asm9s{Xo) 
V2 



1 



A — Ao + ro — iacos^ ro — iacos^ 



P\X) = P\Xo) + 



-P'(Ao) 



^^iasin^ < 
V2 \ 



^'(Ao) 



1 



1 



1 



A — Ao + J'o — iocos^ To — iacosf 
1 



-2aMcos ^s(Ao 



A — Aq + ro + iacos^ ro + iacos^ 
A - Ao + ro 



ro 



(A — Ao + ro)2 + a2 cos2 9 rg + a2 cos2 9 



(63) 
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where sin^ and cos^ refer to the fixed value = Oq. 

For example, on the equatorial plane = 7''/2, Eq. (63) reduces to 

s(A) = s(Ao) , 

52(A) = 52(Ao)-P^(Ao)(A-Ao) + 4sas(Ao)' ^ ^ 



P\X) = P\Xo) + V2aIm[P^(Ao)] ( ^ - -) . (64) 

V A - Ao + ro roj 



In the special case ro = Ao = r+ = M + \/ — and P^(Ao) = 0, i.e., the particle 
starts moving from the horizon r+ , we have A = r and 

s{r) = s{r+) , S2{r) = S,{r+) - P\r+){r - r+) + ^as{r+) - ^ 

P\r) = V2aIm[P^{r+)] - ■ (65) 

If in addition P^{r+) = 0, then we have 

s{r) = s{r+), S2{r) = S2{r+) + ^asir+)^l-^y (66) 

We see that S2 (carrying spin weight —1) is coupled (via the Kerr rotational pa- 
rameter a) to (i.e.,s, with spin weight 0): in a sense, a helicity-rotation coupling 
factor as is necessary to ensure spin-weight conservation on both sides of the equa- 
tion for 5*2. 

The behaviour of the varying component P^ of the total 4-momentum as a 
function of the radial parameter is shown in Fig. 1. P^ is always negative in this 
case, and goes to the asymptotic value — v^aIm[P^(r+)]/r+ for r — > 00. Note that 
in the Schwarzschild limit (a = 0) P^ vanishes identically. 



5.5. The Kasner case 

The Kasner metric^^ representing an anisotropic universe model is 

ds^ = dt^ + t^P'dx'^ + t'^P^dy^ + t'^P^dz^ , (67) 

with pi + P2 + Ps = Pi+ P2+ P3 = 1- Choosing Kasner index values pi = —1/3, 
P2 = Ps = 2/3 (and permutations) one has the Petrov type D case, and introducing 
the NP tetrad 

f-pi fPi 1 

l=—[dt+ t-P'd4 ,n=-^[dt- t-P'd,] , m = ^[t'^^'dy + it'"' 9,] , (68) 

one finds the nonzero spin coefficients 
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while the only surviving Weyl scalar is V2 = — 2/(9f^). The integral curves of I are 



3/2 



1., . .0 t\ 



2/3 



a; = a;o + ^(A-Ao)' + -^(A-Ao) 



t 

y = yo, z = zo , (70) 
and s(A) = s(Ao), P^A) = P\Xo), = P'(Ao), S2 = -P^{Xo){X - Xo) + S2{Xo). 

6. Explicit example in the type N case 

Consider the metric associated with a single gravitational wave^^ in coordinates 
u,v,x,y 

ds^ = -Hdu^ + dudv - dx^ - dy^ , (71) 

where 

H = /+ (u) (x^ + 2/x {u)xy , (72) 

and /+,x represent the two possible polarization states of the wave. The natural 
NP tetrad is 

I = 2dy , n = du+ Tidy , m = A={dx + idy) (73) 

v2 

with the single nonvanishing spin coefficient 

z/ = - [/+ {u) {x + ly) + i/x {u) {x - iy)] (74) 

and Weyl scalar ^4 = — [/+(u) + i/x (u)]. The integral curves of I are 

u = uo, V = 2(A - Ao) + vq, x = xo, y = yo, (75) 

with uo,vo,xo,yo arbitrary integration constants taken at A = Aq. The Mathisson- 
Papapetrou equations for the world line of the multipole reduction chosen within 
the congruence I are 

ds r. dSa dPi „ dP^ 

dA=°' df = -^ ' ^ = «' ^ = °' 

with solution: s(A) = s(Ao), S'2(A) = S'2(Ao) - P^(Ao) (A - Aq), = P^(Ao) and 
p3 = P^{Xo). The case P^ = is trivial. 

7. Null circular orbits in some familiar type D spacetimes 

Instead of choosing Z to be a (geodesic) principal null direction of the background 
spacetime, one can allow it to be the tangent vector to any other null geodesic. Null 
circular geodesies in some of the above type D spacetimes are of special interest; in 
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the black hole case they represent the last possible circular orbits as one approaches 
the horizon. To describe this case, Eqs. (11) reduce to 

dA=0' 



d52 
dA 
dPi 
"dA" 

dA 



TTqIS ■ 



= -ttqP^ - TfoP^ + is {ip2 - V2) 



(77) 



where we have set e = 0, assuming / to be affinely parametrized, and in addition 
tjji = 0, because this is the case in all the explicit examples considered below and 
wc have used the property that all of the nonvanishing spin coefficients (ttq denotes 
the constant value of tt when evaluated along the orbit) and the Weyl scalars are 
constant along the orbit. Then the system can be easily integrated and the solution 
is 

s(A) = s(Ao) , 

P\X) = 2Re{7ro[52(A) - 52(Ao)]} - 2s(Ao) Im[V2(Ao)](A - Aq) + P^(Ao) , 



P3(A) = nois{Xo) - ^ , 
where S'2(A) satisfies the following equation 

-1P0S2 ■ 



d^So 



dA2 



(78) 



(79) 



For example, in the cases in which tjjQ (constant along the orbit) is real and positive, 
the solution for 5*2 is given by 



5*2 (A) = ci cos -y/ VoA + C2 sin ■\/i/'oA 



+ 1 



■ C4e 



(80) 



When instead ^/'o is not purely real, we have to solve a constant coefficient second 
order system of coupled ordinary differential equations for the real and imaginary 
parts of 5'2, which is also trivial to do. We give the essential details to obtain the 
solution for the Kerr, Schwarzschild, Taub-NUT and C metric spacetimes in the 
following subsections. 

7.1. Kerr and Schwarzschild spacetimes 

In the case of the Kerr solution, it is well known that in the equatorial plane 
{9 = 7r/2) there exist two circular null geodesies at radii 

"2 



r-(geo),± = 2M I 



1 + cos 



-arccos 



V m) 



(81) 



which are real solutions of the cubic equation 



6Mr^ + 9MV - AMa^ = . 



(82) 
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To discuss this case, we consider the following NP frame field adapted to generic 
co-rotating and counter-rotating circular null orbits 

2aMrsin6'±\/AS 



l± = d, 



n± = 



m • 



sme[(r^ + a2)Z; + 2a^ Mr sin^ I 
{r^ + a^)j: + 2a^ Mr sin^ 9 



2EA 



2aMrsin6'=F V^S 



sin 6»[(r2 +a?)i: + 2a^Mr sin^ 6*] 



-^(VAdr + ide) . 
\/2E V J 



(83) 



The congruence I contains the null geodesic orbits in the equatorial plane. Along 
such orbits, 9 = tt/2 and r = r(ggo) and the only nonvanishing spin coefficients 
reduce to 

V2{r-M) ^ . VA 



-T = 



l3 = a + 



a = 



V2 (r2 + a2)(r3 - 3Mr^ + ra^ + Ma^) ± aM^/A{'ir'^ + a^) 



4r^ 



and the nonvanishing Wcyl scalars are 

3MA 



V'o 



V'4 



r3(r3 + aP'r + 2a?M) 
3M 



4r.5A 



r'^ +a^± a\/A 



r3 + a^r + 2Ma2 



.n2 



,(84) 



M 



(85) 



where r = r(geo),± is assumed. Thus the spin coefficients and Weyl scalars are all 
constant along the orbit. From Eqs. (f f ) we have 



ds 
dA 



, ^=nis~P-\ 
^ = -n{P' + P'), ^ = ^0^2. (86) 

which can be easily integrated because all the coefficients are constant. 

In the limit of the Schwarzschild spacetime all these results simplify. The NP 
tetrad adapted to the null circular orbits is 



^ ^ Vrjr - 2M) ^ r 
l± = OtT - — :r^—p{ — OA , n± = 



2(r - 2M) 



sm 9 



m 



(87) 



Along the two equatorial null geodesies orbits 9 = tt/2 and r = 3M the only 
nonvanishing spin coefficients are 

V6 



TT = /3 = — T = — 



18M ' 



while the nonzero Weyl scalars 

1 



^0 



'27M2 ' 



^"2 



1 



54M2 ' 



1p4 



(88) 
(89) 



18 Bini D., Cherubini C, Geralico A. and Jantzen R. 



are all constant. From Eqs. (86) we easily get the solutions. 



7.2. Taub-NUT spacetime 

Introduce an NP frame adapted to a generic circular null orbit 
VAtn 



Of. 



(r2 + ^2) sin + 2iy/A^ cos 6 



dr , 



m 



(r^ + e^f sin^ 9 - Atn cos^ i 
V2 



TN 



(r2 +e)shie- 2£v'A^57 cos 9 



VArN^r- + ide 



2(^2+^2)1/2 

Null geodesies exist off the equatorial plane satisfying the conditions 



- 2,Mr^ - 2,fr + = , cos^ 9 



(r2 + £2)2 + 4^2 A, 



TN 



The first equation gives 

r-(geo) =M + 2VM2+£2 

cos 



1 / M 

-arccos = 

3 V\/M2 +£2 



(90) 



(91) 



(92) 



The only nonvanishing spin coefficients (evaluated along the orbits, i.e., using 
Eqs. (91) to simplify general formulas) therefore reduce to 



V2 



a = (3 + 



4(r2 +£2)1/2 

2 (r2+ £2)3/2 
V2 e^/A^ 



r- M 



vAtn 



i cot ( 



4£yA^J7(r2 + £2) cos 61 - [{r"^ + £2)2 + 4£2Atn] sin 9 



4 (r2 + £2)3/2 



(^2 +£2)2 _4^2Atn 



o Atn(»' + £ ) ' sm 6/ 



(93) 



and the nonvanishing Wcyl scalars are 

Atn(^- jM)(£ + ir)3 



^°""^[(r2+£2)4-i6£4A^j,]2l(^ +1^^^' 



4a2 
TN 



+4£A/A^(r2 + £2)[(r2 + £2)2 + 4£2Atn] sin^cosf 



V'2 = 



i £-iM 



M(r2+£2)4_16£4A|N]' , 



2 (r + i£)3 ' 4 A2 j^(r2 + £2)3(£ + irf 

The integration of the equations of motion then follows easily. 



(94) 
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7.3. C metric 

Introduce the following NP frame 

1 



I = -;=[VHdu 
v2 



dr — Ar"^ sin 9d9 



H 



rVCdcj)] 



n 



m 



1 

71 



dr — Ar"^ sin OdO n 
Hdu H ;= r^dcj)] 



dr 



■ r sin 9 



i Ar \ 



Null geodesies exist in the equatorial plane (9 = tt/2) at r = 3M. Evaluated 
along this orbit, the only nonvanishing spin eoefficients are: 

/ M 



while the Weyl scalars reduce to 
tpo = -- ^ 



-a • 



1p2 = 



1 



18M2 ' 54M2 ' 

The equations of motion are easily integrated. 



■^4 = - 



1 



18M2 



(96) 



(97) 



8. Concluding remarks 

We have studied the motion of massless spinning test particles in algebraically 
special vacuum spacetimes according to the Mathisson-Papapetrou model supple- 
mented by Pirani's conditions in the framework of the Newman-Penrose formalism. 
We have considered the case in which the particle follows a null geodesic world line 
belonging to a congruence of a principal null direction of the background space- 
time, a condition which allows complete integration of the equations of motion 
and for this case we have discussed explicit examples for both type D spacetimes 
(Schwarzschild, Kerr, Taub-NUT, C metric, Kasner) and type N spacetimes (grav- 
itational wave), revealing helieity couplings with connection eoefficients as well as 
the curvature of the background. Finally, we have also considered for some familiar 
type D spacetimes the motion of a massless spinning test particle along a circular 
null geodesic, a case which can also be completely solved analytically. 
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